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rR  AC  T 

Hydrostatic  extrusion  is  a promising  new  metal  working  process.  The 
intent  of  this  research  program  is  to  relate  the  process  variables  of  hydro- 
static extrusion  to  the  microstructure  and  properties  of  a wide  variety  of 
extruded  materials  and  to  analyze  this  process  of  metal  forming  by  theoretical 
mechanics.  A hydrostatic  extrusion  press  has  been  designed  and  built  for 
these  experiments.  This  device  is  capable  of  subjecting  a 3/8  inch  diameter 
billet  to  500,000  psi  of  hydrostatic  presure.  Special  tooling  has  been  built 
for  the  mechanical  testing  of  the  experimentally  extruded  materials.  Certain 
speciallized  specimen  materials  preparation  equipment  has  been  set  up,  for 
example  a mechanical  attritor  to  facilitate  the  preparation  of  particulate 
composite  alloys.  The  theoretical  mechanics  analysis  has  been  directed 
toward  the  development  of  a practical  finite  element  analysis  of  the  mechanics 
of  the  hydrostatic  extrusion  process.  Progress  has  been  made  toward  a physically 
and  mathematically  rigorous  formulation  of  the  problem  amenable  to  numerical 
solution  by  means  of  finite  element  computer  methods. 
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Introduction 


Hydrostatic  extrusion  is  the  metal  working  operation  that  consists  in 
forcing  metal  through  a die  aB  in  conventional  extrusion  but  with  the  addition 
of  a superimposed  hydrostatic  pressure.  The  effect  of  the  superimposed 
hydrostatic  pressure  is  to  permit  the  extrusion  of  materials  not  practically 
extrudable  with  the  conventional  process,  Lt  alLows  greater  extrusion  ratios 
to  be  used,  and  it  seems  to  impart  imp: oved  mechanical  properties.  This 
promising  m^tal  forming  technique  has  already  had  some  technological  applica- 
tions, and  its  potential  for  industrial  utilization  is  very  attractive. 

Dispite  the  considerable  technological  promise  of  hydrostatic  extrusion 
there  is  a lack  of  understanding  about  the  basic  engineering  mechanics  of  the 
process,  and  there  is  little  information  about  how  the  variables  of  the  process 
relate  to  the  resulting  mechanical  properties  of  hydrostatically  extruded 
materials.  It  is  the  intent  of  the  present  program  to  make  inroads  on  both 
these  deficiencies. 

Our  work  on  hydrostatic  extrusion  follows  three  main  lines  of  investigation: 

1.  The  systematic  investigation  ol  the  influence  of  the  processing  variables 
(extrusion  ratio,  extrusion  rate,  superimposed  hydrostatic  pressure,  extrusion 
temperature)  and  the  resulting  mechanical  properties  and  microstructure  of  several 
fairly  simple  materials  such  as  pure  Fe,  Ni,  Mg  and  Ti;  and  Cu-30Zr 

2.  Exploratory  attempts  to  improve  the  properties  of  more  complex  materials 
chosen  for  their  technological  utility  such  as  Mg-Li-B  for  high  specific  stiffness, 
Fe-C  for  high  strength,  T-D  nichrotne  for  high  temperature  strength,  Ti-6A1-4V  for 
its  high  strength  to  weight  ratio,  and  NbTi  a metallic  superconductor. 

3.  Analysis  of  the  mechanics  of  hydrostatic  extrusion.  This  will  involve 
both  theoretical  mechanics  and  complementary  experiments  to  test  the  theoretical 
predictions.  The  thrust  of  the  theoretical  mechanics  program  is  toward  the 
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development  of  finite  element  analysis  methods  capable  of  yielding  an  accurate 
numerical  description  of  the  hydrostatic  extrusion  process  taking  into  account 
the  deformation  characteristics  of  the  extruded  materials,  die  geometry,  strain 


rates,  etcetera. 

During  this  first  year  of  work  considerable  progress  has  been  made  toward 
developing  experimental  apparatus  and  techniques,  and  in  overcoming  some  formidable 
analytical  and  computational  problems.  These  are  detailed  in  this  report. 

II*  Experimental  Equipment  and  Techniques 

1.  Hydrostatic  Extrusion  Press 

The  key  piece  of  appartus  for  this  work  is,  of  course,  the  hydrostatic  ex- 
trusion press.  Funds  for  the  purchase  of  the  hyd-ostatic  extrusion  press  were 
provided  by  an  equipment  grant  from  the  National  Science  Foundation  and  by 
Stanford  University.  This  device  has  been  designed  and. built  by  Revere  Copper 
and  Brass  Company.  The  time  r^nuired  lor  the  construction  of  the  extrusion 
press  considerably  exceeded  the  manufac  turers  estimate.  Originally  planned  for 
installation  at  Stanford  University  in  January  of  1975;  it  was  delivered  early 
in  June  1975  and  as  yet  no  experimental  hydrostatic  extrusions  have  been  made 
with  it. 

The  hydrostatic  extrusion  press  was  given  trial  runs  at  the  Rome,  New  York 
pli.nt  of  the  Revere  Copper  and  Brass  Company  where  Mr.  Robert  Whalen  of  Stanford 
observed  and  discussed  all  details  of  the  design  and  operation  of  the  device 
with  Revere  engineers.  Much  of  our  consideration  in  the  initial  operation  of 
the  apparatus  • Lll  be  with  safety,  both  with  respect  to  the  operating  personnel 
and  for  the  equipment  itself.  Safety  is  obviously  a matter  of  great  importance 
in  the  successful  operation  of  ultra  high  pressure  equipment. 

In  most  respects  the  design  and  fabrication  of  the  Revere  Copper  and  Brass 
Company  hydrostatic  machine,  Model  0310,  has  followed  very  closely  the  normal 
pattern  of  prototype  equipment  development.  Major  and  minor  design  alterations 
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and  modlf ications  have  resulted  in  numerous  unpredicted  delays.  It  is  true  that 
one  other  prototype  of  this  basic  design  had  been  built,  but  significant  improve- 
ments have  clearly  established  the  uniqueness  and  superiority  of  the  Stanford 
model.  The  original  machine  should  be  considered  a prototype  for  industrial 
applications  with  lesser  emphasis  on  its  value  as  a basic  research  tool.  Our 
machine  has  been  altered  somewhat  to  provide  more  precise  control  of  those  pro- 
cessing variables  (e.g.  temperature,  extrusion  rate,  bark  pressure)  involved 
with  the  empirical  and  theoretical  predict- Lon  of  material  and  processing  pro- 
perties. Thus  its  utility  is  essentially  for  research  with  indirect  applications 
to  industrial  process  development.  It  is  with  great  enthusiasm  that  we  await  the 
installation  of  the  hydrostatic  extrusion  machine  at  Stanford  scheduled  for 
June  1975. 

The  hydrostatic  extrusion  system  consists  of  three  separate  parts  (Fig.  1) 
the  extruder,  the  hydraulic  control  system,  and  the  control  panel  which  electron- 
ically links  the  two  together.  In  addition  to  the  control  panel,  the  extruder 
and  hydraulic  system  are  interconnected  by  a system  of  microswitches  mounted  on 
the  extruder.  These  switches,  which  will  be  discussed  in  more  detail  later, 
establish  the  position  of  the  piston  for  extrusion  and  behave  as  limit  switches 
as  a precaution  against  over-extruding. 

The  hydraulic  system  was  custom  manufactured  by  Airline  Hydraulic  Corporation 
to  provide  the  necessary  oil  flow  and  pressure  as  regulated  by  the  specific  function 
controls  to  be  set  by  the  operator.  Since  the  extruder  and  hydraulic  system  work 
together,  it  is  best  to  combine  the  discussion  of  their  operation. 

From  the  standpoint  of  design,  the  extruder  is  the  more  interesting.  The 
extrusion  cylinder  and  mounting  stand  occupy  a volume  of  approximately  sixty-two 
inches  long  by  sixty-four  inches  high  by  forty-eight  inches  wide.  The  extrusion 
cylinder  which  is  about  fifty-two  inches  long  is  dixtdr.d  into  two  parts  with 
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Figure  1.  Hydrostatic  Extrusion  system 


distinct  diameters.  The  smaller  diameter  section  of  the  cylinder  contains  the 
extrusion  piston  and  die  stem  guides.  On  top  of  the  cylinder  are  located  the 
previously  mentioned  microswitches  (set  Fi>;.  2 and  drawing  3a).  Also  visible 
in  Fig.  2 are  the  hydraulic  ports  mounted  on  the  face  and  wall  of  the  small 
diameter  cylinder.  The  attached  larger  diameter  port  of  the  cylinder  contains 
the  extrusion  load  chamber  and  all  the  accompanying  equipment  necessary  for  the 
containment  of  the  high  pressures  which  develop  during  extrusion  (see  Fig.  2 
and  drawing  3b). 

The  floating  die  jacket  design  for  the  containment  of  the  high  hydrostatic 
pressures  is  quite  interesting.  It  was  chosen  over  other  designs  for  several 
reasons.  For  one,  the  cylinder  fatigue  life  (estimated  to  be  106  cycles)  is  far 
superior  than  with  other  systems  designed  to  hold  comparable  pressures.  Secondly, 
this  design  is  more  tolerant  of  dimensional  changes  which  must  naturally  occur 
under  such  loading  conditions.  When  changes  do  occur  the  machine  will  either 
continue  to  operate  normally  or  require  only  minor  adjustments  of  the  die  size. 
Finally,  loading  and  reloading  for  either  direct  or  indirect  extrusion  is  rapid, 
permitting  many  extrusions  in  a day. 

Although  some  of  the  ideas  and  principles  involved  are  ingenious  the 
operation  of  the  extrusion  press  is  quite  simple.  To  extrude  in  the  indirect 
mode  the  following  steps  are  taken.  The  reverse  button  is  pushed  on  the  control 
panel  to  fully  retract  the  piston,  die  stem,  and  die  stem  guides.  The  micro- 
switches are  adjusted  for  the  proper  extrusion  stroke  length.  Next  the  die, 
billet,  oil  and  high  pressure  seal  placed  Into  the  loading  tube  shown  in  Fig.  4. 
The  loading  tube  is  used  to  transfer  the  assembled  die,  billet,  oil,  and  seal 
intact  into  the  liner  of  the  die  Jacket  in  the  high  pressure  extrusion  chamber. 
Once  the  transfer  has  taken  place,  a tungsten  carbide  anvil  is  secured  against 
the  seal  to  prevent  backward  motion  due  to  the  hydrostatic  pressure,  (see  drawing 


Figure  3a.  Cross  sectional  drawing  of  extrusion  piston,  die  stem 
and  die  stem  guides. 


Figure  4.  Proper  loading  tube  assembly  for  indirect  extrusion.  In 
the  order  shown,  die,  billet,  oil,  and  high  pressure  seal  are  loaded 
into  the  tube  at  the  bottom  of  the  photo. 


Figure  5.  Cut-away  view  of  loaded  hydrostatic  extrusion  chamber  - 
indirect  extrusion  mode.  The  material  has  been  transferred  intact 
from  the  loading  tube  (Figure  3). 
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3b,  and  Figs.  5 and  6.  Note  the  direction  in  which  the  material  extrudes.)  The 
machine  is  now  ready  for  the  extrusion  cycle.  A push  of  the  "pressure  intensifier" 
button  directs  a system  hydraulic  pressure  of  3,000  psi  to  i.he  intensifier  piston 
(located  approximately  half  way  up  the  load  frame,  see  Fig.  2).  The  pressure 
intensifier  boosts  the  pressure  to  60,000  psi  and  delivers  the  high  pressure 
through  stainless  steel  high  pressure  tubing  to  the  load  chamber  cylinder,  pres- 
surizing the  external  surface  of  the  die  jacket.  The  uniform  compressive  force 
or.  the  die  jacket  outer  surface  performs  two  functions.  First,  it  provides  the 
necessary  counterbalancing  compression  force  on  the  liner  to  contain  the  high 
internal  pressures  developed  during  extrusion.  Secondly,  it  contracts  the  liner 
sufficiently  to  form  a metal  to  metal  seal  with  the  die,  a condition  necessary  for 
the  prevention  of  oil  leakage  from  the  load  chamber.  By  activating  the  extrusion 
function  at  the  control  panel,  the  piston  and  die  stem  begin  to  move  at  a pre- 
selected rate  toward  the  load  chamber  until  the  hollow  die  stem  contacts  the  die. 

At  this  point  the  extrusion  pressure  will  build  until  the  material  begins  to 
extrude.  Provided  the  necessary  pressure  to  extrude  is  not  greater  than  the  max- 
imum system  pressure,  the  piston  will  always  move  at  the  selected  rate.  If  the 
pressure  to  extrude  is  too  great,  a static  condition  will  prevail,  and  no  extrusion 
will  take  place.  Other  values  for  the  processing  variables  must  then  be  chosen. 
Since  all  the  load  is  carried  by  the  die  stem,  the  die  stem  guides  are  a necessary 
precaution  against  buckling  of  the  die  stem.  When  the  piston  has  travelled  a 
certain  pre-set  length,  a microswitch  is  activated  which  automatically  stops  the 
extrusion.  This  is  a precaution  against  over-extruding,  a potentially  costly 
error  which  will  be  discussed  later.  Removing  the  materials  from  the  load  chamber 
after  extrusion  is  a simple  procedure.  First,  remove  the  breach  plug  securing  the 
anvil  against  the  load  chamber  seal.  Then  remove  the  anvil.  Now,  by  using  an 
over-ride  button  on  the  control  panel  move  the  piston  in  the  forward  direction 
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until  the  seal,  remaining  billet  material,  and  die  are  pushed  out  of  the  load 
chamber  by  the  die  stem.  Protection  from  moving  the  piston  too  far  forward  and 
jamming  it  against  the  high  pressure  cylinder  is  provided  by  another  microswitch 
which  8 huts  the  system  off  when  tripped. 

The  same  procedure  is  followed  for  extrusion  in  the  direct  mode,  except  that 
the  order  and  direction  in  which  the  material  (seal,  material,  die)  are  loaded  is 
reversed.  See  Figs.  7 and  8 for  the  positioning  in  the  loading  tube  and  die  jacket 
liner,  respectively.  Note  also  that  the  extrusion  direction  is  opposite  to  in- 
direct extrusion. 

The  machine  is  interchangeable  for  either  three-eights  inch  or  three-quarter 
inch  billets.  To  change  the  system,  the  die  jacket  assembly,  die  stem  and  die 
stem  guides  must  be  removed  and  replaced  by  the  other  set.  Such  an  exchange  con- 
sumes about  three-quarters  of  a day.  The  machine  is  designed  to  operate  at  maximum 
system  pressure  (3,000  psi)  with  safety  when  rigged  for  a three-quarter  inch  billet. 
A hydrostatic  pressure  of  250,000  psi  is  developed  in  the  load  chamber  when  the 
piston  cylinder  pressure  is  3» 000  psi.  One  must  exercise  caution  when  using 
the  three-eights  inch  billet  system  since  hydrostatic  pressures  of  1,000,000  psi 
could  be  developed  with  a piston  cylinder  pressure  of  3,000  psi,  while  both  the 
die  stem  and  die  jacket  were  designed  to  contain  pressures  only  up  to  500,000  psi. 
For  this  reason,  the  maximum  permissible  extrusion  pressure  developed  by  the 
hydraulic  system  is  adjustable;  and  it  is  only  a ~atter  of  a few  calculations  to 
determine  the  appropriate  maximum  allowable  extrusion  pressure  before  critical 
stresses  build  up  in  the  die  stem  and  the  die  jacket. 

The  hydrostatic  machine  is  a precision  engineered  design  with  most  tolerances 
±.001  inches.  With  proper  care  and  respect  no  significant  wear  should  occur. 
Fortunately,  most  moving  parts  are  wel]  lubricated.  The  cit  of  die  size  to  die 
jacket  liner  size  must  be  constantly  checked  to  ensure  metal  to  metal  contact 
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during  extrusion.  Since  the  die  and  liner  do  not  make  well-] ubricated  contact, 
some  wear  is  expected.  The  liner  should  periodically  be  checked.  However,  it 
is  constructed  of  polished  tool  steel  and  usually  will  be  harder  than  the  chosen 
die  material.  Therefore,  its  lifetime  should  be  somewhere  around  250  to  1,000 
cycles.  Its  replacement  is  of  minor  consequence.  If  a die  is  found  to  be  in 
good  condition  but  undersized,  the  diameter  may  be  built  up  easily  by  electroplating 
with  copper. 

As  with  most  high  pressure  equipment  potentially  hazardous  and  costly  situa- 
tions may  arise.  All  precautions  will  be  taken  from  the  start  to  ensure  that  the 
operator  is  well  protected  from  any  possible  machine  failure.  All  high  pressure 
cylinders  and  tubing  will  be  covered  by  heavy  metal  shielding.  The  main  pressure 
cylinder  should  not  present  a major  problem  since  the  material  was  chosen  for  its 
ability  to  absorb  energy  through  slow  crack  propagation.  The  volume  of  oil  com- 
pressed at  60,000  psi  around  the  die  jacket  is  too  small  to  impart  much  energy  to 
the  system  during  a failure.  In  any  case,  these  areas  will  be  protected  by  shields. 
More  importantly,  the  machine  will  be  oriented  so  that  both  ends  of  the  extrusion 
cylinders  face  barriers.  It  is  quite  possible  under  certain  conditions  for  high 
pressure  oil,  die  stems,  or  material  projectiles  to  be  ejected  at  dangerous 
velocities  from  either  end. 

Care  must  be  exercised  to  prevent  two  things  that  could  go  wrong  which  would 
cause  no  physical  harm  but  would  require  costly  replacement  or  repair.  (1)  The 
intensifier  pressure  could  drop  below  60,000  psi  resulting  in  a split  liner  and 
expanded  die  jacket  assembly,  and  (2)  the  die  stem  could  break. 

The  die  stem  could  be  broken  by  simply  miscalculating  and  exceeding  the 
compressive  yield  stress  of  the  die  stem  during  extrusion  with  the  three-eights 
inch  system.  Failure  of  the  die  stem  might  also  result  from  a sudden  loss  of 
hydrostatic  pressure  in  the  die  jacket  load  chamber.  This  might  occur  from 
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overestimating  the  extrusion  stroke.  At  the  beginning  of  the  stroke,  the  pressure 
builds  until  the  extrusion  is  initiated;  steady  state  extrusion  follows.  But 
instead  of  triggering  a microswitch  to  halt  the  forward  motion  of  the  piston  before 
the  billet  is  completely  extruded,  the  piston  continues  forward  and  extrudes  the 
entire  billet  creating  a pressure  drop  to  atmospheric  pressure  in  the  chamber. 

Such  almost  instantaneous  unloading  on  the  die  stem  could  cause  critical  tensile 
stresses  to  develop  and  the  die  stem  to  fracture.  This  condition  may  also  arise 
if  oil  loss  occurs  in  the  load  chamber.  Extrusion  of  oil  instead  of  material  due 
to  improper  billet  to  die  seal,  or  cracked  dies  can  result  in  such  a failure. 

However,  in  most  situations,  loss  of  oil  can  only  stall  the  machine  at  maximum 
extrusion  piston  pressure  with  no  damage  resulting. 

Failure  or  loss  of  intensifier  pressure  is  potentially  a very  serious  problem. 
Machine  damage  (although  reparable)  is  almost  certain  to  occur  if  the  pressure  on 
the  outer  surface  of  the  die  jacket  drops  below  some  critical  value  during  extrusion. 
This  critical  pressure  varies  according,  to  the  applied  extrusion  pressure  and  whether 
the  three-eights  inch  or  three-quarter  inch  system  is  in  use.  As  a safety  measure, 

the  die  jacket  is  always  pressurized  to  the  maximum  allowable,  i.e.  60,000  psi. 

Failure  could  result  in  several  ways.  Perhaps  the  most  serious  failure  mode  would 
be  when  the  pressure  seals  actually  prevent  the  high  pressure  oil  from  entering  the 
clearance  between  the  die  jacket  and  outer  large  diameter  cylinder  (see  Fig.  9a, b 
for  seal  design).  The  pressure  in  the  line  and  on  the  gage  would  read  the  equiva- 
lent 60,000  psi,  but  the  actual  pressure  within  the  area  of  the  die  jacket  would  be 

near  atmospheric  pressure  (see  Fig.  10a, b).  The  other  potential  mode  of  failure  is 
through  leaks  in  the  seals  or  in  the  check  valves  in  the  hydraulic  line. 

Because  of  the  seriousness  of  the  problem,  we  will  mount  strain  gages  externally 
on  the  large  cylinder  over  the  area  of  the  die  jacket.  The  strain  gages  will  be 
connected  by  a relay  network  to  the  control  panel.  The  machine  will  be  programmed 
not  to  operate  if  either  pressure  does  not  develop  around  the  die  jacket  or  if 
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Figure  10b.  Possible  seal  failure.  Due  to  the  absence  of  a high  compensating 
compressive  stress,  tensile  stresses  now  develop  and  .likely  lead  to  load 
chamber  liner  failure. 


pressure  suddenly  drops  around  the  die  jacket. 

Much  of  the  initial  work  done  with  the  machine  will  involve  calibration  of 
the  strain  gages,  flow  meters  and  determination  of  die  friction  and  oil  compressi- 
bi1 ity.  Presently , the  applied  load  to  the  die  is  measured  as  strain  on  a strain 
gage  resulting  from  the  load  transmitted  to  it  from  the  die  stem.  This  is  calibra- 
ted with  t.ie  pressure  in  the  extrusion  piston  cylinder  to  give  a numerical  value 
of  the  applied  load  and  therefore  hydrostatic  pressure.  Although  this  technique 
accurately  measures  the  load  on  the  die  stem,  one  must  be  careful  in  relating  it 
to  the  hydrostatic  pressure  in  the  load  chamber  or  the  pressure  at  the  extrusion 
die.  In  the  indirect  mode  of  extrusion  the  most  accurate  results  of  hydrostatic 
pressure  at  the  die  are  possible.  Since  the  die  stem  and  the  die  are  in  contact, 
only  the  wall  friction  between  die  and  liner  contributes  to  the  difference  between 
internal  hydrostatic  pressure  and  load  on  the  dLe  stem.  Through  experiments  it  will 
be  possible  to  gain  an  approximate  value  for  die  friction  under  various  condition. 

The  direct  mode  of  extrusion  presents  considerable  problems.  The  die  stem  is  no 
logger  close  to  the  extrusion  die  and  under  high  pressures  the  oil  tends  to  solidify, 
thus  gaining  the  capacity  to  carry  a load.  This  result  is  a pressure  drop  from  the 
point  of  application  of  the  load  (i.e.  the  die  stem),  to  the  actual  point  of  ex- 
trusion. Such  readings  at  the  die  stem  obviously  involve  unpredictable  errors.  It 
may  be  possible  to  mount  strain  gages  on  the  die  for  the  direct  extrusion  mode  to 
get  more  reliable  results. 

Oil  compressibility  in  the  load  chamber  is  something  which  should  be  known  in 
order  to  accurately  position  the  microsswitches  for  proper  extrusion  stroke  length. 
Knowledge  of  the  compressibility  will  become  expecially  important  when  using  the 
back  pressure  system.  Here  the  available  length  of  the  outlet  chamber  into  which 
the  material  will  be  extruded  will  be  related  to  fluid  compressibility.  This  length 
will,  or  course,  vary  with  pre-set  back  pressure  reading. 
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While  it  is  expected  that  che  back  pressure  system  will  be  an  extremely 
useful  tool  in  the  future,  it  is  arriving  at  Stanford  as  untested  equipment. 

Such  a system  has  never  been  used  before  with  this  machine.  Even  thouph  the 
back  pressure  system  is  a simple  design  it  may  require  considerable  trial  and 
error  to  make  it  function  properly. 

2.  Testing  Equipment 

Mechanical  property  measurements  will  be  an  important  part  of  this  program. 

It  is  planned  that  both  tension  compression  and  fatigue  tests  will  be  performed 
on  hydrostatically  extruded  experimental  materials.  For  that  purpose  a fully 
instrumented,  20,000  pound  capacity,  hydraulic  mechanical  testing  system  (MTS 
model  No.  901.29,  Research  Inc.)  has  been  set  up  adjacent  to  the  location  of  the 
hydrostatic  extrusion  press.  This  mechanical  test  equipment  is  fully  operational. 

In  addition  to  room  temperature  properties,  we  are  also  interested  in  mech- 
anical properties  at  elevated  temperatures.  Furthermore,  we  are  restricted  to  very 
small  compression  samples  due  to  the  limited  size  of  the  extruded  materials.  To 
meet  these  needs,  a furnace-cooling  water  system  has  been  fitted  onto  the  MTS  frame 
in  such  a way  that  we  can  compress  very  small  samples  (as  small  as  0.15  inches  in 
diameter  by  0. 200  inches  long)  in  air  up  to  1000°C.  This  furnace  attachment  was 
not  simple  to  set-up  because  of  two  design  problems.  The  first  problem  was  the 
non  uniform  deformation  of  the  small  compression  samples  caused  by  even  the 
slightest  misalignment  of  the  compression  columns.  This  problem  was  eventually 
solved  (or  at  least  minimized)  by  decreasing  the  length  of  the  compression  columns 
to  6 inches.  This  in  turn  required  that  the  furnace  be  placed  on  rollers  and 
allowed  to  slide  along  a vertical  track.  The  second  problem  was  the  conduction  of 
heat  out  of  the  furnace  and  into  the  load  cell  and  hydraulic  fluid  on  the  MTS  frame 
by  the  compression  columns.  This  problem  was  eliminated  by  installing  alumina 
columns.  The  furnace  MTS  system  is  now  fully  operational.  A photograph  of  the 
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system  is  included  as  Figure  11. 

3.  Materials  Preparation  Equipment 

Included  in  our  planned  hydrostatic  extrusion  program  is  the  study  of 

hydrostatically  extruded  particulate  composites.  Such  composites  are  oftentimes 

difficult  to  extrude  by  ordinary  methods  and  should  be  ideally  suited  for  testing 

the  feasibility  of  hydrostatic  extrusion.  Earlier  work  by  us  in  the  preparation 

of  particulate  composites  involved  the  use  of  a special  powder  raetallurgy-mech- 

(1  2) 

anical  comminution  process'1  ’ developed  at  Stanford.  This  work  hcs  been 
successful  in  obtaining  unusual  mechanical  properties  in  large  volume  fraction 
particulate  composites  based  on  model  systems  (zinc  and  cadmium  base  metals  con- 
talnlng  alumina,  tungsten  and  boron)  . Recently  we  have  concentrated  on 
magnesium  and  copper  base  particulate  composites.  As  discussed  in  the  first 
semi-annual  report ^ these  systems  have  posed  some  new  exDerimental  difficulties. 
1)  The  compaction  press  previously  used  could  only  deliver  a maximum  of  30,000  psi 
pressure  on  a 1 1/2  inch  bore  die,  which  is  inadequate  for  hot  pressing  of  these 
newer  materials.  2)  A large  size  difference  between  the  matrix  powders  (30-40 
microns)  and  the  second  phase  powders  (about  1 micron)  prevented  uniform  distribu- 


tion of  the  hard  particles  in  the  matrix  material.  Copper  and  magnesium  readily 
oxidize  at  the  elevated  temperatures  of  hot  pressing  which  inhibits  interpnrticle 
diffusion  and  results  in  a poorly  bonded  final  product  (hence  cracks  during  ex- 
trusion) . 

The  first  problem  has  been  solved  by  the  use  of  a 120,000  pound  capacity 
pneumatic  compaction  press  provided  by  the  Civil  Engineering  Department  at  Stanford. 
This  virtually  doubles  our  maximum  compaction  pressure. 

The  second  problem  has  been  solved  by  the  acquisition  of  a mechanical  attritor 
from  Union  Process  Co.,  of  Akron,  Ohio.  The  attritor  not  only  mixes  powders  very 
well  but  also  can  grind  powders  in  a protective  atmosphere  to  submicron  size. 
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Figure  11.  Photograph  of  the  MTS  model  No.  901.29,  Research  Inc. 
with  the  corresponding  high  temperature  mechanical  testing  system 
attached  to  the  load  frame.  Not  shown  is  the  furnace  temperature 
controller  and  the  water  cooling  system  for  the  load  cell. 


Figure  12.  Photograph  of  the  attritor  acquired  from  Union  Process  Co. 
as  installed  in  our  Powder  Metallurgy  Laboratory. 
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Largo  particle  matrix  powders  can  thus  be  reduced  to  a size  comensurate  to  the 
second  phase  particles.  Our  attritor,  only  recently  delivered  (May  9,  1975),  has 
already  been  installed  and  used  to  grind  cast  iron  powders.  The  first  test  runs 
have  only  been  partially  successful  - the  powders  are  still  coarse  even  after 
many  hours  of  grinding.  This  can  be  attributed  to  the  large  size  of  the  charged 
cast  iron  particles  (0.010  inches  thick  and  0.25  inches  long).  Better  results  are 
expected  as  the  influence  of  the  process  conditions  on  the  effectiveness  of  the 
mechanical  comminution  are  learned  through  experience.  A photograph  of  the 
installed  attritor  is  included  as  Figure  12. 

The  third  problem  has  been  solved  by  the  design  of  an  inert  atmosphere  hot- 
pressing  furnace,  which  is  shown  schematically  in  Figure  13.  This  furnace  allows 
for  the  hot  pressing  of  powders  at  hi*h  stresses  (up  to  60,000  psi)  and  temperatures 
(up  to  900°C)  in  a protective  (or  reducing)  atmosphere.  The  compaction  chamber 
is  floating,  pressure  being  applied  simultaneously  to  both  the  top  and  bottome 
oi  the  powder  charge,  thereby  minimizing  density  gradients  along  the  length  of 
the  compacts.  The  device  is  presently  under  construction. 
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Figure  13.  Schematic  diagram  of  the  inert  atmosphere  hot-pressing  furnace 
currently  under  construction.  This  furnace  is  designed  to  handle  up  to 
900°C  enviornments  and  compaction  stresses  up  to  60,000  psi. 


— ■ 


1 


I 


III.  De ve  1 opra ent  of  a Flnlte-Elomei it  P rogra m for  Met al  Forming 


It  was  pointed  out  in  the  First  Semi-Annual  Technical  Report,' ^ that 
in  spite  of  the  large  strains  which  occur  In  metal  forming  processes, 

It  is  necessary  to  carry  out  an  elastic-plastic  analysis  in  order  to  be 
able  to  anticipate  the  limits  on  satisfactory  operation.  Even  though 
the  plastic  strains  are  very  large  compared  with  elastic  strains,  the 
latter  play  the  major  role  in  determining  the  stress  distribution,  and 
this  governs  the  development  of  defects,  through,  for  example,  the  initi- 
ation of  internal  fractures  and  instabilities.  Elastic  strains  are  of 
the  order  (yield  stress) / (e last ic  modulus),  which  is  about  10'3  for  most 
technologically  important  metals.  Thus  computational  programs  must  be 
accurate  to  this  order  of  strain  if  stresses  are  to  be  accurately  predicted 
throughout  the  work-piece.  In  regions  where  large  plastic  flow  is  cur- 
rently taking  place,  the  stress  is  effectively  determined  by  the  plasticity 
laws,  but  in  the  final  product,  for  example,  plastic  flow  has  ceased  and 
elastic  strains  determine  the  distribution  of  residual  stresses  which 
are  important  from  a technological  standpoint.  Thus,  while  it  might 
appear  from  a purely  formal  computational  standpoint  that  since  strains 
of  the  order  unity  occur,  strains  of  magnitude  10'3  could  be  neglected, 
this  is  not  the  case  for  metal  forming  analysis  aimed  at  assessing  satis- 
factory operational  limits  for  process  variables. 

The  need  to  adopt  elastic -plastic  theory  has  a major  influence  on 
the  mathematical  structure  of  the  analysis.  This  aspect  of  the  p-oblera 
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has  received  concentrated  attention  of  the  theoretical  group  during  the 
first  year  of  the  project.  Sotre  of  the  considerations  which  arise  in 
this  connection  are  discussed  in  subsequent  sections  of  this  report. 
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1.  Rigid-Plastic  and  Elastic-Plastic  Theories 

For  simplicity  of  i-he  initial  discussion,  consider  ideal  plasticity 
obeying  the  Mises  yield  condition.  The  stress  at  a material  particle 
is  given  by  the  Cauchy  or  "true"  stress  tensor  a , with  six  independ- 
ent components  because  it  is  symmetric:  the  direct  stresses  , 022 

and  and  the  shear  stresses  ct^  » and  . It  is  convenient 

to  write  the  plasticity  laws  in  terms  of  the  stress  deviator,  4^  _ , 
which  is  the  stress  tensor  less  the  average  hydrostatic  component: 


*ij  = aU  3°kk6i  j 


(1) 


where  6 is  the  Kronecker  delta: 


6ij  “ 1 * 1 


* 0 , i 4 j 


(2) 


and  the  summation  convention  for  the  repeated  suffix  k is  utilized. 

The  deviator  expresses  plastic  flow  characteristics  since  it  in  effect 
represents  the  shear  type  stress  components.  Plasticity  occurs  when  the 
yield  stress  is  reached  at  the  Mises  limit: 


where  Y is  the  yield  stress  in  tension.  The  equations  determining  the 
plastic  deformation  are  flow  or  incremental  in  nature,  giving  the  rate 
of  plastic  strain  cj  , or  equivalently  the  increment  in  plastic  strain 
in  a time  increment  dt 


(4) 


For  ideal  or  perfect  plasticity,  flow  occur.1  with  the  stress  remaining 
at  the  yield  limit  (3),  and  the  plastic  strain-rate  is  given  by 


(5) 


1 


where  X in  a free  parameter,  since  indefinite  flow  can  occur  at  constant 


stress  with  ideal  plasticity  and  thus  no  work-hardening.  Since  6^  . is 


a deviator  in  (5),  so  also  must  gj  be,  hence  the  rate  of  volume  strain 


is  zero: 


•Ji-« 


(6) 


These  relations  govern  the  plastic  components  in  elastic-plastic  theory, 
to  which  must  be  added  the  elastic  components  of  strain-rate: 
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(7) 


(8) 


(9) 


where  G Is  the  shear  modulus. 


Equations  (5)  and  (7)  Illustrate  the  different  structures  of  rigid- 


plastic  and  elastic-plastic  theories.  For  rigid-plastic  theory,  e 


•P 


ij 


is  equal  to  the  total  strain-rate  « . For  a velocity  field  v(x,t) 


the  subscript  curls  denote  vectors  - where  x are  position  coordinates 
in  a body  in  the  current  (or  deformed)  configuration,  the  strain-rate 
is  given  by 


l/*Vi 


3v 


2 V 3x , 


3x, 


(10) 


For  a given  velocity  field,  (5)  in  conjunction  with  (3)  determines  the 


stre?s  deviators,  4 ^ , for  rigid -plastic  response  to  stress,  whereas 


for  elastic-plastic  analysis,  (8)  yields  only  a differential  equation 


for  the  determination  of  4^  , or  equivalently  it  yields  4^  if  4^ 


is  known.  This  difference  is  well  illustrated  by  the  graphical  construction 
for  the  case  of  principal  directions  of  stress  and  strain  fixed  in  the 
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body,  known  as  the  tt  diagram/8^  Figures  14  and  15. 

The  tt  diagram  is  a projection  of  the  principal  stress  space  onto 
the  plane 

al  + °2  h a t " 0 (11) 

that  is  a projection  parallel  to  the  line  subtending  equal  angles  with 
the  three  stress  axes.  Strain  can  be  plotted  on  the  same  diagram,  and 
the  scale  factor  2G  , where  G is  the  elastic  shear  modulus,  produces 
coincidence  of  the  stress  and  strain  points  for  purely  elastic  deformation 
Fig.  la  shows  the  rigid-plastic  case  for  which  the  plastic  strain-rate 

* P 

e is  identical  with  the  total  strain  rate.  The  strain-rate  vector 
DE  is  parallel  to  the  stress  vector  OA  . Thus  given  the  strain-rate 
DE  , the  stress  point  A on  the  yield  circle  is  uniquely  determined. 

In  the  elastic-plastic  case,  Fig.  lb,  DE  again  represents  the  total 
strain-rate  vector,  now  broken  down  into  a plastic  component  FE  and 
an  elastic  component  DF  . The  latter,  including  the  factor  2G  , is 
equal  to  the  stress-rate  vector  AB  . In  this  case,  the  total  strain- 
rate  vector  DE  does  not  determine  the  stress  OA  , but  given  OA  , it 
then  determines  the  plastic  component  of  the  strain-rate  FE  parallel 
to  OA  , and  hence  the  elastic  component  DF  . 

(9) 

Prandtl  gives  a geometrical  or  mechanical  analogy  for  determining 
the  stress  variation  associated  with  a prescribed  strain  path.  For  ideal 
plasticity,  if  a head  is  dragged  behind  the  strain  point  on  a thread  of 


length  Y , the  thread,  as  the  bead  follows  a curve  of  pursuit,  gives  the 
stress  vector  according  to  (8).  The  curve  of  pursuit  in  effect  solves 
the  differential  equation  (8). 

The  two  types  of  response  to  prescribed  jtrain-rate  just  illustrated 
both  analytically  and  graphical ly,and  corresponding  to  rigid-plastic  and 
elastic-plastic  behavior,  represent  quite  different  ways  in  which  the 
history  of  straining  influences  the  current  stress.  For  rigid-plastic 
material,  the  rate  of  strain,  or  more  precisely  the  velocity  strain  (10), 
determines  the  stress  deviator  uniquely.  The  stress  tensor  itself  is 
determined  to  within  an  arbitrary  hydrostatic  pressure,  which  is  associ- 
ated with  the  incompressible  characteristic  of  plastic  flow.  Thus  no 
influence  of  the  prior  history  of  strain  appears  in  the  equations  which 
determine  the  stress.  In  contrast,  for  the  elastic-plastic  case,  the 
current  stress  is  given  by  the  solution  of  the  differential  equation  (8) 
integrated  over  the  history  of  deformation  from  prescribed  initial  condi- 
tions. This  circumstance  determines  the  structure  of  the  solution  process 
for  stress  and  deformation  distributions  in  metal  forming.  Consider  the 
solution  to  have  been  obtained  up  to  the  time  t . The  stress  and  defor- 
mation distributions  are  then  known  at  that  instant.  The  equilibrium 
equations,  constitutive  relations  and  velocity  strain-rate  relations  can 
be  most  conveniently  written  in  terms  of  stress  rates,  for  as  in  the 
case  of  Equation  (8)  these  are  determined  uniquely.  Equivalently,  for 
a sufficiently  small  time  step  dt  , stress,  strain  and  displacement 
increments  are  determined.  The  resulting  new  stress  and  deformation 
distributions  constitute  the  "initial  conditions"  for  the  next  time  step. 
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An  analogous  situation  arises  in  the  case  of  more  realistic  work- 
hardening laws  of  plasticity,  which  again  require  integration  through 
time  for  the  determination  of  stress  distributions.  A generalized  tt 
diagram  construction  has  also  been  dwe loped  for  this  case  ^8»10). 

It  is  interesting  to  observe  that  the  extent  to  which  strain  history 
affects  the  current  stress  has  a bearing  on  the  convection  of  stress 
magnitudes  through  motion  of  the  material  elements.  In  the  rigid-plastic 
case  this  would  be  zero  since  the  current  velocity  strain  distribution 
determine?  stress  quite  uninfluenced  by  its  prior  values.  In  the  actual 
evaluation  for  elastic-plastic  materials,  this  aspect  appears  in  the 
appropriate  choice  of  convected  derivative  lor  substitution  in  the  con- 
stitutive relation. 


The  need  to  evaluate  the  stress  and  deformation  distributions  in 
terms  of  stress  rates  and  velocities  in  each  time  increment  raises  the 
question  of  the  appropriate  rates  for  expressing  the  equilibrium  equations 
and  the  constitutive  relation.  The  former  can  clearly  be  most  simply 
treated  in  terms  of  the  partial  tire  derivative,  for  this  commutes  with 
the  space  derivative  in  the  equilil  rium  equation: 


?fij 

9 


= 0 


(12) 


according  to  the  relation: 


JL^LL 

?)  t V 9xJ 


= -l/lZii 

a x l a t 


(13) 


However  the  stress -rate  term  in  the  constitutive  equation,  . , must 
fol  ow  the  material  particle,  and  Rice  ^hns  pointed  out  that  however 
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small  the  time  increments,  substitutiot  of  cr  for  the  partial  derivative 
with  respect  to  t in  (13)  would  usually  involve  major  error. 

If,  for  illustration,  is  talen  to  be  the  material  derivative, 

then  (13)  takes  the  form 


A.  (*_I±± 

dM  9 1 


A./-  daii 

N ' v»  ^ 


(14) 


For  a time  increment  dt  , we  write 


d<Tij  ■ V'  • <*«!,  - , d„  - y dt 


(15) 


and  (14),  multiplied  by  dt  , takes  the  forn 


dT  (dCTlj) 


(16) 


In  order  of  magnitude,  d <j  ^ E t d ; , wheru  E*  is  the  appropriate 
tangent  modulus,  (du^)  - de  , so  i hat  the  first  term  dominates  if 


E » .r 


(17) 


L icrti 
In  this  case  it  is  permissible  to  repl  ice  -r— 4 by  a. . and  the  influence 

o t i j 

of  convection  on  the  material  stress  rite  is  negligible.  For  loading  in 
the  plastic  range  of  many  metals,  the  tangent  modulus  is  of  the  order  of 
the  stress,  so  that  the  terms  in  (16)  ire  of  the  same  order.  Thus  the 
circumstance  which  applies  in  classical  elasticity,  for  which  the  elastic 
modulus  is  large  compared  with  the  striss,  that  the  difference  in  using 
deformed  or  undeforraed  coordinates  is  legligible,  does  not  usually  apply 
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for  plastic  loading.  It  is  clear  from  (15)  and  (16)  that  this  circum- 
stance is  independent  of  the  size  of  the  time  increment  dt  . It  is 
perhaps  worth  pointing  out  that  the  additional  terms  comuonly  added  to 
generate  objective  stress-rates  are  of  the  same  order  as  the  second  term 
in  (16)  so  that  the  same  conclusion  results.  Moreover,  when  steady  state 
operation  of  a metal  forming  process  has  been  achieved,  3a^/9t  ■ 0 , 
and  the  two  terms  in  the  second  expression  in  (14)  are  exactly  equal 
in  magnitude,  whatever  the  characteristics  of  the  material,  and  must 
clearly  both  be  retained. 

2.  Formulation  of  Finite-Element  Relations 

A discussion  of  the  generation  of  finite  element  relations  for  metal 

(12) 

forming  problems  has  been  given  by  l.ee  and  Kobayashi  and  forms  the 
basis  for  the  present  discussion,  home  of  the  pitfalls  which  must  be 
avoided  in  applying  this  approach  to  elastic -plastic  analysis  are  empha- 
sized. 

In  accordance  with  the  structure  of  elastic-plastic  theory  discussed 

in  the  previsou  section,  which  showed  that  stress  and  strain  rates  (or 

equivalently  increments  of  stress  and  strain  in  the  time  At)  are  sought 

on  the  basis  of  knowledge  of  the  current  stiess  and  deformation,  a var- 

(13) 

iational  principle  in  the  stress  and  strain  rates  is  utilized  : 

6$  = 0 (18) 


where 
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Superscript  T denotes  the  transpose  ol  a matrix  or  a vector.  The  stress 
tensor  a is  written  as  a vector  o with  nine  components,  as  also  is 

the  strain  tensor,  so  that  the  first  integral  gives  the  rate  of  work  carried 
out  by  the  stress-rates.  The  second  integral  expresses  the  rate  of  external 
work  expended  by  the  traction  rate  vector  T acting  over  the  surface  S . 

The  body  being  analysed  is  divided  into  M finite  elements  considered 
to  be  interconnected  at  a finite  number  of  nodal  points  P . The  func- 
tional (19)  evaluated  for  the  m1"*1  element  is  denoted  by  <t^m\  so  that 


The  discretization  of  the  variational  problem  for  each  element  is  accom- 
plished by  approximating  the  functional  by  a function  cj/m^  of 

the  m1"*1  element  nodal  point  values.  This  is  achieved  by  replacing  the 
v distribution  with  an  approximating  velocity  distribution  in  each  element, 
continuity  between  elements  being  maintained.  The  functional  <t>  , (19), 
is  thus  approximated  by  a function  of  the  P nodal-pcint  velocity  vari- 


ables v 
~P 


* 2 ^i*  v**  V = S <p(m)  (v£m)) 


22. 


Tc- 


Thus  the  Initial  variational  problem  is  reduced  to  the  determination  of 

Lite  vector  v in  the  P-dimensiona  1 vector  space  of  vectors  which  min- 
-P 

imize  cp(v  ) . This  vector  combined  with  the  elemental  velocity  distribu- 
-P 

tions  provide  the  approximation  to  the  solution  for  the  actual  velocity 
distribution  v(x,t)  . 

The  approximate  velocity  variation  in  the  m^  element  is  given  by 


— (m)  , v r.. , (m)  (m) 

vv  (x , t ) = [N]  ( x)  v (t) 


(22) 


where  v and  [ N ] are  of  the  m^  element  nodal  point  velocity 

components  and  the  mt'1  distribution  matrix  respectively.  From  (22)  the 
approximate  strain-rate  vector  can  he  expressed  in  terms  of  the  nodal 
velocities  by  the  matrix  operator  [ B ] ; 


(tn) 


|B]<">  v'"1 


(23) 


Both  in  the  elastic  and  plastic  regions,  the  vectors  of  rates  of  stress  and 
strain  are  related  by 


(24) 


(n) 

where,  for  elastic  response  the  matrix  [D]  has  constant  elements,  and 
for  plastic  behavior  it  depends  on  the  current  stress  and  on  the  loading 


history.  Substituting  (22),  (23)  and  (24)  in  (19)  for  the  m 
gives 


th 


element 


;[D](m)fBl(m)v(ra))T[B](n,)  v(m)  dV 


1 


<[Nl<m,v,m>)'r  ds 


3 


Evaluation  of  the  volume  and  surface  integrals  for  a single  element 
results  in  the  discrete  representation  of  the  functional  for  an  element 


cp (m)  = i v("0T[K](m)  v(m)  . V(*)T  ‘(m) 


where  [K]  is  the  element  stiffness  matrix  and  R m'  is  the  nodal 

force-rate  vector  equivalent  to  the  external  traction  rate  T distri- 
buted over  the  surface  S(m)  . The  global  stiffness  matrix  [K|  and  the 
equivalent  nodal  force-rate  vector  R ire  the  sums  of  the  corresponding 
element  values  as  called  for  in  (21).  The  condition  for  cp  to  assume  a 
stationary  value  yields  a system  of  simultaneous  linear  algebraic  equations 
which  determine  the  node  velocities  v^  . From  these,  stress-rates  can  be 
obtained  from  (24)  and  (23)  . 

When  mixed  boundary  conditions  are  considered  for  which  velocity 
vectors  over  part  of  the  surface  S(;  are  prescribed,  the  unknown  nodal 
velocity  vectors  are  given  by 


R * [K]  v + | K]  , v 
-a  aa  .a  Jab  _b 


where  R is  the  known  nodal  force-rate  vector  and  v and  v are 
~a  ,a  _b 

respectively  the  unknown  and  prescribed  nodal  velocity  vectors.  The 


24. 


matrices  [1<1  and  [K]  are  siibmatrlces  of  the  global  stiffness  matrix 
na  ab 


[K] , according  to 


[K]  v = 


ba  i 1 Jbb  ( lb 


! 


Note  that  in  this  finite-element  formulation  the  stress  rate  appearing  in 

the  functional  (19)  was  obtained  from  the  material  constitutive  relation 

(24)  through  the  relation  (23)  based  on  nodal  point  velocities.  Thus 

rates  of  stress  and  strain  following  specific  particles  are  evaluated, 

called  material  derivatives  (convect ed  derivatives  in  the  fluid  mechanics 

/ 1 3 \ 

literature).  But  (19)  was  established'  from  a virtual  work  hypothesis 
which  assumed  that  0^  satisfied  equilibrium  equations  of  the  form  (12). 
We  saw  in  the  last  section  that  this  i'-  not  valid  for  plastic  flow  since 
the  partial  time  derivative  30^/St  is  not  approximated  by  the  material 
derivative,  and  it  is  the  former  which  satisfies  the  equilibrium  equation. 
The  situation  becomes  particularly  anomalous  in  the  case  of  steady-state 
problems,  since  then  the  partial  derivative  of  the  stress  is  zero,  and 
substitution  of  the  material  derivative  for  it  is  totally  suspect. 

As  discussed  by  McMeeking  and  Kice^^  , Hill^^  has  attempted  to 
avoid  this  difficulty  by  working  with  stress  as  force  per  unit  reference 
area  (Piola-Kirchhof f stress)  instead  of  the  "true"  or  Cauchy  stress 
utilized  so  far  in  this  report,  which  is  defined  per  unit  current  area. 
McMeeking  and  Rice  consider  an  increment  of  deformation  from  a current 
state  which  for  the  time  interval  At  is  taken  to  be  the  reference  state. 
For  plasticj ty  laws  it  is  not  appropriate  to  use  a permanent  reference 
state  to  define  the  stress  and  strain  tensors,  since  plastic  flow  is  an 


incremental  or  flow  type  phenomenon  which  can  most  naturally  be  expressed 
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in  the  current  geometry  as  for  a fluid^  . However,  taking  the  current 


states  sequentially  as  a series  of  icference  states  for  a sequence  of  de- 


-mation  increments,  flow  type  response  can  be  accommodated  with  the 


advantage  of  a referential  theory  rather  than  a spacial  theory  (Euler 
coordinates).  In  a referential  type  theory  (Lagrange  coordinates)  the 


equations  of  motion  are  written  in  terms  of  the  reference  coordinates,  so 
that  partial  differentiation  with  respect  to  time  corresponds  to  a fixed 


reference  or  material  point,  and  so  is  simultaneously  a material  rate. 


Even  when  an  appropriate  reference  formulation  is  adopted,  care  must 
be  exercised  in  selecting  appropriate  tensor  rate  definitions,  for  consti- 


tutive relations  must  be  objective,  and  this  is  most  easily  achieved  if 
all  variables  occurring  in  them  are  objective.  Objectivity  is  the  require 


ment  that  a constitutive  relation  f..r  a material  must  be  such  that  if  n 


time  dependent  rigid  body  rotation  is  superposed  on  the  body,  the  stress 
determined  by  the  constitutive  equation  for  the  new  motion  must  simply  be 
the  original  current  stress  rotated  by  the  current  value  of  the  superposed 
rotation.  In  this  definition  it  is  assumed  that  any  change  in  inertia 
forces  is  equilibrated  by  body  forces  so  that  they  do  not  complicate  the 


situation. 


Some  variables  which  would  appear  to  be  plausible  candidates  for 


constitutive  relations  are  found  not  to  be  objective  and  are  therefore 


unsuitable.  For  example  the  materi il  time  derivative  of  stress 


• do , . d(  , , 

o = — * 1 + v 

ij  3t  Vk  dx. 

K 


is  not  objective  (see,  for  example,  Malvern  p.  402).  But  the 


co-rotational  stress-rate  or  Jaumann  rate,  which  is  equal  to  the  material 
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rate  less  the  stress  rate  associated  with  th<*  rotation,  is  objective  [10]. 
Thus  an  appropriate  rate  is  needed  (or  substitution  In  the  elastic-plastic 
relation  (24),  and  the  influence  of  this  choice  on  the  form  of  the  vari- 
ational principle  (L9)  must  be  investigated. 

We  are  looking  into  these  questions  to  cevelop  the  most  effective 
form  of  finite-element  analysis  on  which  to  base  our  attack  on  the  steady- 
state  extrusion  problem.  We  have  already  formulated  possible  approaches 
to  the  satisfaction  of  steaay-state  boundary  conditions  at  fixed  special 
positions  as  the  work  material  moves  through  them.  We  expect  to  run 
corresponding  numerical  programs  dui ing  the  summer. 
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